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I Abstract. This paper generalizes the static and spherically sym- 

metric traversable wormhole geometry to a rotating axially sym- 
C , metric one with a time-dependent angular velocity by means of an 

exact solution. It was found that the violation of the weak energy 
condition, although unavoidable, is considerably less severe than in 
the static spherically symmetric case. The radial tidal constraint 
is more easily met due to the rotation. Similar improvements are 
^ ■ seen in one of the lateral tidal constraints. The magnitude of the 

oo ■ angular velocity may have little effect on the weak energy condi- 

' tion violation for an axially symmetric wormhole. For a spherically 

symmetric one, however, the violation becomes less severe with in- 

■ creasing angular velocity. The time rate of change of the angular 
I velocity, on the other hand, was found to have no effect at all. 

■ Finally, the angular velocity must depend only on the radial coor- 
^ , dinatc, confirming an earlier result. 
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1. Introduction 



^ . It was recognized by Flamm (jy in 1916 that our universe may not 

be simply connected: there may exist handles or tunnels, now called 
wormholes, in the spacetime topology linking widely separated regions 
of our universe or even connecting us with different universes alto- 
gether. That such wormholes may be traversable by humanoid travel- 
ers was first conjectured by Morris and Thorne [J^, thereby suggesting 
that interstellar travel and even time travel may some day be possible. 
For a detailed discussion see the book by Visser [Sj. 

Morris-Thorne (MT) wormholes are static and spherically symmetric 
and connect asymptotically flat spacetimes, here assumed to be isomet- 
ric. Adopting units in which c = G = 1, the metric for this wormhole 
is given by 

(1) ds^ = -e^^Mrft^ + + r'idO^ + sin2^t/02); 

1 — b[r) jr 
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$(r) is called the redshift function and b{r) the shape function. The 
shape function describes the spatial shape of the wormhole when viewed, 
for example, in an embedding diagram, described below. 

To hold such a wormhole open, violations of certain energy conditions 
proved to be unavoidable. More precisely, all known forms of matter 
obey the weak energy condition (WEC) Tapn^ii^ > for all timelike 
vectors and, by continuity, all null vectors (Friedman Matter that 
violates this condition is called exotic by Morris and Thorne. 

Various attempts have been made to generalize the MT wormhole 
by giving up spherical symmetry [3j or by including time-dependence. 
A particular interesting example of the latter is the inclusion of a de 
Sitter scale factor multiplying the spatial part of the metric. The goal 
was to study the possibility of enlarging a wormhole pulled out of the 
spacetime foam to macroscopic size (Roman 5J). A similar scale factor 
was used by Kim |6| . Yet another possibility is the use of a conformal 
factor n{t) [3 Em GDI: 

(2) ds' = m [-e^^^'^dt^ + e'^'^'-^dr^ + r\de^ + sin^^ d<p')\ . 

For a metric with time-dependent functions $ and A, see Ref. All 
these studies include a discussion of the WEC violation. Traversabihty 
conditions are investigated in T] and TT]. 

In this paper we generalize the MT wormhole in another direction 
by assuming the wormhole to be rotating, not necessarily at a constant 
rate, and by dropping the assumption of spherical symmetry. Instead, 
the wormhole is assumed to be axially symmetric, i.e., symmetric with 
respect to the axis of rotation. Stationary axially symmetric wormholes 
are discussed by Teo jT^j. It is shown that the WEC is indeed violated 
but that a traveler would not necessarily come into contact with any 
of the exotic matter. These wormholes may also have an ergoregion, 
where a particle cannot remain stationary with respect to spatial in- 
finity. This is an extreme example of the well-known dragging effect in 
general relativity. 

The main purpose of this paper is to discuss both the energy viola- 
tion and the traversabihty conditions by first finding an exact solution. 
Possible restrictions on the metric coefficients recently proposed by 
Perez Bergliaffa and Hibberd ^S] are discussed in Section [7| 

Proposals to search for naturally occurring wormholes, if they exist, 
go back at least to 1995 jT3]. For a summary of these findings see 
Ref. While definite conclusions are still lacking, the possible exis- 
tence of wormholes or the existence of negative mass cannot be ruled 
out. 
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2. The metric 

The study of spacetimes that are both stationary and axially sym- 
metric has a long history El E]- A spacetime is stationary if 
it possesses a time-hke Kilhng vector field rj°- = [d/dt)"" generating 
invariant time translations. Axially symmetric is formally defined as 
possessing a space-like Killing vector field = {d/d(f)Y generating 
invariant rotations with respect to 0. 

Suitable metrics for stationary axially symmetric fields are discussed 
by Islam [T^. We will adopt the metric suggested by Teo |12^, since it 
appears to be best-suited for obtaining an exact solution: 

(3) ds'^ = -N^dt^ + e^c/r^ + r^iT^ [dO'^ + sin2^(rf0 - u dtf] , 

where N, /i, K, and u are all functions of r and 6; uj is the angular 
velocity d(j)/dt. More precisely, 

_d(j) _ d(j)/dT _ 
dt dt /dr ' 

referred to in Ref. [201 as the "angular velocity relative to the asymp- 
totic rest frame." The connection between the metric due to a bounded 
rotating source and the mass and angular momentum of the source is 
discussed in Ref. jTH]. 

To make our solution as general as possible, we will assume that 
UJ = uj{r,6,t) is time-dependent, so that the wormhole can no longer 
be called stationary. The reason for proposing this model is a practical 
one: if an advanced civilization were to succeed in constructing such a 
wormhole, it is likely to be aspherical and the rate of rotation likely to 
be varied. Accordingly, we will write our line element as follows: 

(4) = -e^^^^'^^dt^ 

+ e^'^^'-'^^c/r^ + [K{r, e)f [dO^ - sin^e{d(P - u{r, 6, t)dty] . 

Here K{r, 9) is a positive dimensionless function of r such that Kr 
determines the proper radial distance at (r, 9) in the usual manner. 
In other words, 2Ti{Kr)sm9 is the proper circumference of the circle 
through (r, 9). 

So far nothing has been said about the shape function. Recall that 
the wormhole geometry may be conveniently described by means of 
an embedding diagram in three-dimensional Euclidean space at a fixed 
moment in time and for a fixed value of 9, the equatorial slice 9 = 
7r/2 |2l 1^. The resulting surface of revolution has the parametric 
form 

f{r,(f)) = {r cos (f),r sin (p, z{r, 9 1)), 
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where z = z{r, 6) is some function of r and 6 and 6 is momentarily 
held fixed at 6i. As usual, we think of the surface as connecting two 
asymptotically flat universes. The radial coordinate decreases from 
+00 in the "upper" universe to a minimum value r = tq at the throat, 
and then increases again to +oo in the "lower" universe. 

Since the embedding surface must have a vertical tangent at the 
throat for any value of 6', we require that 

dz 

lim — = +00, 

while \imr-^oo dz/dr = 0, the meaning of asymptotic flatness. Return- 
ing to the line element (j3]), we further assume that for any fixed 6, 
fi{r,9) has a vertical asymptote at r = tq: \imr^ro+ fJ'{r,9) = +oo. 
Also, fJ>{r,9) is a twice different iable function of r and 9 (as is \{r,9)) 
and is a strictly decreasing function of r with limr^oo f^{f, = 0. 
These requirements are met hy z = z{r, 9) for any fixed 9 such that 

dr 

(for the upper universe). Furthermore, d'^z/dr'^ < near the throat 
(since dfj.{r,9i)/dr < 0), as required by the "flaring out" condition in 
Ref. [2]. The shape function is now defined by 

1 _ fc(ne) • 

r 

It follows that 

Finally, at the throat itself, b must be independent of 9. It is readily 
checked that db/d9 = at r = ro. 

3. The solution 
Let us write the line element (jH) in slightly more compact form: 

(5) ds^ = -e^^dt" + e^^c/r' + if [d9^ + sin^9{d(P - u dtf] . 

To make the analysis tractable, we choose an orthonormal basis {cq,} 
which is dual to the following 1-form basis: 

(6) 0° = dt, 9^ = dr, 9^ = Kr d9 
and 



(7) 9"^ = Kr sin9{d(t)-ujdt). 
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(See also Ref. |23-) a result, 

(8) dt = e-^ dr = e-^ e\ dO = -^6^, 

Kr 

and 

(9) d(p = —^e^ + uj e-"" e\ 

Kr sm 

Furthermore, 

ds" = -{de^f + {de^f + {de^ + {de^f. 

Since the orthonormal basis is itself rotating, some information re- 
garding uj is going to be lost. We will briefly return to this topic at 
the end of the last section. (In particular, it will be shown that the 
magnitude of u must be restricted.) To do so, we need the compo- 
nents of the fundamental metric tensor {gap}-, as well as {(7"^}, in the 
(t, r, 6*, 0)-coordinate system: 

(10) gu = -e^^ + irV^cu^ sin^^, gt^ = -K'^r^uj sin^O, 



grr = e^^, gee = K^r"^, g,j,^ = K'^r'^ sw?9, 

and 

(11) (^" = -e-2\ (7*^ = -cue-2\ g" = e-^^', 

The last component, g^'^, bears a striking similarity to d(f) in Eq. Q. 

To obtain the curvature 2-forms and the components of the Riemann 
curvature tensor, we use the method of differential forms (Ref. j22j). 
To that end we calculate the following exterior derivatives in terms of 

d9' = ^e^^ 9^A90 + ^^e'A 9\ d9' = ^^9'A9\ 
or Kr o9 Kr o9 

d9^= (^e-^^ + L^e-^^) 9^A9\ 
\r K or J 

and 

ill ) /// ) 

d9^ = -Kr—e-^e-''sm9 9^ A9^ - — e'^sin^ 0^ ^ qo 
or 39 

The connection 1-forms ^ have the symmetry 



1, 2, 3), and = -uj\ {t, j = 1, 2, 3, ^ ^ j] 
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and are related to the basis 6^ by 

The solution of this system is found to be 

2 or or 
^ Kr 09 2 89 



}:Kr^e-'e-'^sm9 9' + -%-\in9 9', 
2 dr 2 39 



^ ^ Kr 89 V r ^ ~^ K dr^ 



^1 = \Kr^e-^e~^'sm9 9"^ - ( -e'^ + -^e''^ 
2 dr \r K dr 

uj\ = -_e-\in0 0° - — cot0 



2 09 \Kr K^r 09 

The curvature 2-forms ^2*^- are calculated directly from the Cartan 
structural equations 

The results are given in the Appendix. Since the components of the 
Riemann curvature tensor can be read off directly using the formula 

n« ^ * /J*" A an- 

" j - -2^mnj ^ /\ ^ ) 

there is no need to list them explicitly. As an example, suppose we let 
m = and n = 1 in the equation 

Then 

= -2-^011° A 6^^ - -^101° A 6*° 



= -i?oii° 9^ A9^ = Ron 0^ A 

Thus -Rgii^ ~ ^(1)0) ill the Appendix. (As in Ref. [22], we omit the 
hats whenever numerical indices are used.) 
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4. WEC VIOLATION 

As with any traversable wormhole, we would expect a violation of 
the weak energy condition near the throat: Tj^/x"/i'^ > for all null 
vectors. As in Morris and Thorne P| and Roman we use a radial 
outgoing null vector /x" = (/i*, /i*", 0, 0) = (1, 1, 0, 0). In our orthonormal 
frame we expect to have the usual stress-energy components T^f, Tff, 
Tgg, T^^, as well as T^^, which represents the rotation of the matter 
distribution |12ij . 

From G^^ = R^^ - jRg^^ and Rat = Racb^ we have Goo + Gu = 

RoO + Rll = -Roii° + -^022° + -^033° " -^011° ~ -^122^ ~ -^133^- ^ short 

calculation yields 



(12) 



2 

= -e 
r 

1 

+ 



2^t 



dX ^ 9/i 

dr dr 



2 dK 
K dr 



-2^1 



+ 



1 



dX 

m ~ 

/dX 

+ 



cot I 

de 

2 d^K 



^2^2 



dr 



(9^2 902 

4 dK 



dX 
dr 
d^X 



Kr dr 



2/t 



-2/i 



duj\ 

de) ' 



1 I ) e~^hiYi'e 



(We have used the Einstein field equations = SvrT, 



-a/3- 



To put this rather long expression in perspective, consider the static 
spherical case jTTj, where only the first term survives: 



(13) 



1 

8^ 



- I e 

r 



2n 



dX ^ dfi 
dr dr 



Referring to Sect. |2l recall that the throat corresponds to the value 
r = Tq. So if /i is a smooth function of r and given that limr^r.o+ f^{f) = 
+00, it follows that \imr^ro+ d^/dr = — oo. In Ref. [llj, A(r) = —n/r, 
K > 0, so that dX/dr = n/r"^, making p — t negative near the throat. 
That a violation of the WEC cannot be avoided regardless of the choice 
of A(r) can be seen geometrically. For if limr^ro+ dX/dr is positive and 
finite, then the sum on the right side of Eq. (|13|) is negative near the 
throat. This might be avoided if limr^ro+ d^/dr = +00. But then 



lim^ 



A(r) 



-00 and e^^*^^-* 0, which yields an event horizon. 
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The main goal in this section is to show that for a rotating axially 
symmetric wormhole the WEC violation is, in principle, much less se- 
vere than for the static spherically symmetric case by a suitable choice 
of A and /i. The functions were chosen primarily for convenience, keep- 
ing the analysis simple and accommodating the next section at the 
same time. 

Taking k, A, and B to be constants, let 

3 



(14) Xir,e) 



k 






r 


h 





TT 



0<0 <-, k>0, 



(15) fi{r,9) 



0]+B 



o<e< 



TT 

2' 



k,B>0, 



where A is large enough to keep the expression inside the brackets 
positive, and 

ke r/TT 
r — ro L V 2 

where e > is a small constant. On the interval [7r/2, vr) each function 
is defined to be the "mirror image," i.e., A and fi are symmetric about 
6 = 71/2. So the discussion may be confined to the interval (0,7r/2]. 
The partial derivatives are listed next for easy reference: 



dX 

dr 
dX 

m 



k 



2k 



1 
3 

71 

2 

7C 

2 ' 



TT 



+ A 



e 



(9/i 

dr 

m 



ke 



[r - roj 
ke 



e]+B 



r 



0. 



On the right side of Eq. (jl2j) . the first two terms are similar to those 
in Eq. (fTSj) . The first term is therefore negative. The fourth term is 
stricly positive for 6 ^ 7r/2, while the third term is positive near the 
throat. In the fifth term the expression 

2 / o \ 2 



2^2 



(r - ro)2 

is similar to the first term but the quadratic factor reduces the size 
of the second term, thereby making the fifth term less harmful. 

To study the effect of K{r,6), we adopt a function similar to one 
suggested by Teo [12], 

(4a sin6')2 



Kir, 



1 



Now the second and sixth terms are positive, as well. The seventh 
term is close to zero near the throat, thanks to the factor e~^'^, and so 
is completely overshadowed by the third and fourth terms. 
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The last term is more of a problem, being strictly negative. However, 
for Teo's choice, uj = 2a/r^, where a is the total angular momentum of 
the wormhole, the last term is zero. Requiring u to be independent of 6 
may be unavoidable, a conclusion also reached by Khatsymovsky 
who states that for a macroscopic wormhole to exist, the angular ve- 
locity must be independent of 9. 

With the last term eliminated, we see that a drastic reduction in the 
energy condition violation is indeed possible, at least in principle. 

Remark: We will see in the next section that the absolute value of 
the first term exceeds that of the second term. Otherwise the WEC 
violation would appear to have been eliminated completely. 



5. Traversibility conditions 



Another area in which improvements over the static spherically sym- 
metric case are possible is in the study of tidal constraints. In partic- 
ular, for an infalling radial observer the components of the Riemann 
curvature tensor are found relative to the following orthonormal basis 
(from the usual Lorentz transformations): 

(16) 

Cq, = 7e£ T 7 (^) Cf , ep = Tief + 7 (^) e£, e^, = e^, eg, = e^. 

A traveler should not experience any tidal forces larger than those on 
Earth. As outlined in Ref. j2], the radial tidal constraint is given by 



l-Ri'o'i'o'l — 



9(i 



1 



c2x2m (108 m)2' 
assuming an observer 2 m tall. We have 



(17) l-Ri/Q'i'o'l — 1-^1 



ftft I 



-2fi 



dr"^ 



dX dfx 

dr dr 



+ 



dr 



+ 



1 dXdfi 



In the static spherically symmetric case only the first term survives. 
The second term, which involves the angular velocity, reduces the size 
of \Rriri\- the first term. 
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is positive near the throat because the positive middle term inside the 
brackets contains the factor (r — tq)^ in the denominator. For the same 
reason the first term is larger than the absolute value of the second 
(near the throat). Since the second term is negative, the net result, so 
far, is a reduction in the size of |-R^m|. 

Concerning the last term, we need to remember that our wormhole is 
not spherically symmetric. The tidal forces experienced may therefore 
depend on the direction of approach. So we must ask the traveler 
to approach the throat in the equatorial plane 9 = 7r/2, and here 
d\/d9 = 0. 

To study the first of the lateral tidal constraints, |-R2'6'2'6'I — (10^m)~^ 
we have from Eq. (fT^ . 



I 2 I I 2 f^\'^ 



As usual, this is a constraint on the velocity of the traveler. It is 
assumed in Ref. that the spaceship decelerates until it comes to rest 
at the throat. So only the first term needs to be examined: 



:i8) \R 



■eiet\ ~ 

IdX 
r dr 



R, 



■022 



1 



+ 



1 dKdX 
K dr dr 



1 dK d\ 



K^r^ dO dO 



3 (dujV 

l\de) 



2A 



sin^^ 



In the static spherically symmetric case only the positive first term 
survives. Once again requiring that the traveler approach the throat 
in the equatorial plane 9 = 7r/2, the next three terms are zero. The 
last term is also zero due to the earlier requirement du/dO = 0. The 
fourth term is negative, and for our choice of K{r, 9) 



1 dKdX 



K dr dr 



< 



IdX 
r dr 



We therefore have a reduction in the size of iRaiat 



met I • 
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For the remaining lateral tidal constraint we need to examine 



(19) 



R 



'<l>t<i>t 



033 



IdX 

r Or 

2 



1 dKdX 
K dr Or 



-2fi 



1 dKdX 1 fduj\ 



- 



K^r^ de 06 A\de J 



e 



-2A 



sin^^ 



The first two terms appeared in the other lateral constraint; the result 
is a reduction in the size of 



R 



■<j>t<j)t 



. Unfortunately, the next term is pos- 
itive. Although small due to the factor e~^^, the overall result is hard 
to quantify and may actually be less favorable than the corresponding 
static case. (It helps that the next two terms are zero for 9 = 7r/2, 
while the last term vanishes due to the requirement duo/ 89 = 0.) 



6. The effect of rotation 

In proposing a model with a time-dependent a;, i.e., an angular ac- 
celeration or deceleration, it was hoped that the findings in the last 
two sections would be strengthened. As it turns out, however, none of 
the derivatives with respect to time occurring in the curvature 2-forms 
found their way into the earlier calculations. 

This failure suggests that the effect of rotation be studied from a 
different perspective. Since we are using a rotating basis, some infor- 
mation regarding uj was lost: no w's appear in any of the curvature 
2-forms, although the derivatives do. So it may be useful to examine 
the weak energy condition violation relative to the (t, r, 9, 0)-coordinate 
system. 

The expression for Tu + T^.^, calculated by the traditional method 
using the fundamental metric tensor (Equations (fTUj) and (fTT|)). contain 
the old terms in Eq. (|12j). but not, of course, in the orthonormal frame. 
The new terms all contain uo: 



SniTtt + Trr) = old terms + R, 
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where 



-new 



2 




i?new obviously vanishes if uj does. The last four terms vanish if we 
assume, as before, that du/dO = 0. With the functions used earlier, 
the first term is negative and the second positive. Unless dK/dr is 
very small, this result is, once again, hard to quantify. 

The situation is rather different if we return to the assumption of 
spherical symmetry, while retaining uj, assumed to be positive. Then 
K = 1 and A and fi are independent of 6. In the [t, r, 6, 0)-coordinate 
system, using the outgoing null vector (1, 1, 0, 0), we have 



Since uj > 0, the last term is positive. So if uo is large, the WEC 
violation is much reduced. 

While the resulting reduction in the WEC violation is a welcome 
surprise, some words of caution are in order. The last term in Eq. 

suggests that if u is large enough, the WEC violation can be 
eliminated altogether. But as explained by Teo jT^j, for a rapidly 
rotating wormhole it may not be possible to use a radially outgoing 
null vector since the gu component of the fundamental metric tensor 
may no longer be negative, as can be seen from the line element, Eq. 



It was pointed out in a recent paper by Perez Bergliaffa and Hibberd 
[TB] that the metric (@)) used in this paper may require further restric- 
tions. In particular, it is shown that a wormhole of the type studied 
by Teo jT2] cannot be generated by a perfect fluid or by a fluid with 
anisotropic stresses. In the first case the condition 



(20) 




7. Additional considerations 



(21) 



G 
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is violated and in the second case, 



(22) 



Goo + > 2 Gi 



03- 



These conditions seem to be met if the metric (j3)) includes the func- 
tions A and fi used in this paper, that is. Equations (fT^ and (fT^ . 
respectively. 

For example, to check Condition (j2T| . a short calculation yields 



Each term contains e"'^, so that G12 ~ near the throat. 

Even more favorable is the outcome of the check on Condition ()22|1 : 
five of the terms in G03 contain the factor duo/ 86, which is equal to 
zero. The remaining terms all contain the factor e~^'^, so that Goa ~ 
near the throat. The left side. Goo + ^33, contains only two terms that 
are neither strictly positive, nor zero, nor contain the factor e~^'^. But 
these two terms are strongly overshadowed by several terms that are 
strictly positive. So Condition (j22|) is easily met in the vicinity of the 
throat. 

Since the vicinity of the throat is the only region that really matters, 
we can construct (in the usual manner) a solution with a suitable radial 
cut-off of the stress-energy tensor. Our solution can then be joined, at 
least in principle, to an external solution also satisfying the desired 
conditions. (For a discussion of the required junction conditions, see 



It is indeed surprising that A and yU, which were chosen for entirely 
different reasons, are sufficient for satisfying Conditions (j2H) and (j22|) . 
thereby overcoming the objections raised to the wormhole in Teo jl2j . 
Unfortunately, Ref ^2] discusses other conditions, some of which are 
not so easily checked. So it may still be necessary to "incorporate more 
realistic (in the astrophysical sense) features, for example heat flux" 
(Ref PS]). Future investigations of these issues could be extended to 
include stability conditions, especially for rotating wormholes, since 
such conditions may also require "more realistic features." 




1 d^X 
Kr drdO 



1 dKdK 
K^r dr 09 



Ref. [2ni.) 
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8. Conclusion 



In this paper the MT wormhole solution was generahzed to worm- 
holes which are both rotating and axially symmetric, i.e., symmetric 
with respect to the axis of rotation. It was concluded that the unavoid- 
able violation of the weak energy condition is less severe than in the 
spherically symmetric case. The radial tidal constraint is more easily 
met due to the rotation. An improvement was also found in one of the 
lateral tidal constraints. Making the angular velocity uj time-dependent 
does not help since none of the time-derivatives in the Appendix ap- 
peared in the calculations. Furthermore, uj must be independent of 
9, in agreement with Ref. [2H]- Finally, the magnitude of the angular 
velocity may have little effect on the WEC violation for an axially sym- 
metric wormhole. In the spherically symmetric case, however, a rapid 
rotation will result in a reduction in the WEC violation, as long as uj 
is not excessively large. 



= 0) e^^e^ + a(2, o) e'^^e^ + a(3, o) e'^^e^ + a(3, i) e'^ a 

+ ^(3,2) 0^ A^^ 



Appendix 



This appendix lists the curvature 2-forms VL"" ■. 



where 




A(2,0) 




e 



1 dKdX 
K'^r dr 89 



e 
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A{3, 1) = -^K^e-'e-^'^sme- ^Kr^e-^-^^^inO 
^ ' 2 dr 2 dr^ 

1 duj ^ 2/, • /I 1 9uj da ^ o,, . „ 

2 or or 2 or or 

I doj dn . . ^ 3 dK duj n,, . „ 

1 dKduj _x -u ■ n 1 -X -a ■ n 

-e e ^sm9 ^ „„ e e ^smb*. 



K de dr 2 drdO 

Q% = B{1, 0) A ^° + B{2, 0) ^2 A ^0 + fi(3, 0) 6^ A 

+ B(3, 1) e^Ae^ + B{3, 2) ^3 A e\ 

where 0) = A(2, 0), 

Rfom--— -^/^ 1 1 /^Ay 1 dKdX 



rdr K^r^de"^ K'^r^ \de J K^r^ 86 86 



2 86 8r 2 8r 86 2 8r 

1 8K8uj ^ 1 8K 8uj _ x _„ . ^ 
e e ''sm^ K7r~^^ ^ '^sm^, 



2K 8r 86 2K 86 8r 

3 8uj ^ r^du} _2„ . „ 1 8K 8uj _2,, . „ 
-e cos^ K—e e ^sm6 r-:——e e ^sm6 



2Kr86 2 8r 2 8r 8r 



1 8K8uj . ^ 
-e smt^. 



XV 86 86 

Q°3 = C(l, 0) A ^° + C(2, 0) ^2 A e° + C(3, 0) A e° + C(2, 1) A 
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where 

C(1,0) = —Kr^ ( ^] e-^V^sin^, 



2 dt\dr 
1 d fdu 
'2di K'de 



C(2,0) = --^f^^) e-^^sin^, 



r or K or or 4 \or J 



I d\ „ 1 1 /(9cu^ ^ 



^ , , 1 (9u; (9A A _„ . ^ 1 (9u; a _,, . ^ 
1 9a; , „ . ^ 1 A /, 

1 dK duj _A _„ . „ 1 OK duj X -„ ■ n 



2K Or de 2K 06 dr 

Q\ = D{3, 0) e^Ae° + D{2, 1) e^Ae\ 

where £'(3,0) - -C{2,1), 

2 (9fr 9„ 1 (9/i 9„ 1 1 fda^^ 



ar rdr K^r^ 39^ K^r^ \ dO 

X is: 9r dr K^r^ 09 09 ' 

n\ = E{1, 0) 9^A9°+E{2, 0) 9^A9^+E{3, 1) ^^Ae^+E(3, 2) 9^A9^, 
where 0) = -A{3, 1), E(2, 0) = -B{3, 1), 

m 1) = A^e-'' - l^e-^ + l^e-- - l^^e-'^ 
' i^'r (9r r Or K Or"^ K Or Or 

, 1 „ 1 OKO^Ji „ 1 a^TT 

^ ' ' Kr^ 09 Kh- Or 09 K^r Or09 

^ ^ J^^.^^ '^^ ^-2X^-1X^:^^20 1 OK OK 

4 Or 09 K^r Or 09 

= F(l, 0) 9^A9^+F{2, 0) A^°+F(3, 1) 9^A9^+F{3, 2) 9^A9^, 
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where 0) = -A{3, 2), F(2, 0) = -B{3, 2), F(3, 1) = E{3, 2), 

1 1 S^ir 1 air , 1 fduj\^ 2A . 2., 



K^r^ \de J K^r"^ 09^ K^r^ 89 A \ d9 J 
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